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Introduction to Seminar on Materials Science Issue 


957Q0020A Moscow IZVESTIYA ROSSIYSKOY AKADEMII NAUK MEKHANIKA 
ZHIDKOSTI I GAZA in Russian No 5, Sep-Oct 94 pp 3-4 


[Introduction for issue: "From the Editor"] 


[FBIS text] Mechanics under conditions of weightlessness and 
space-based materials science, which have been growing vigorously 
in the last 20 years, have now reached a point at which many 
methodological approaches -- as well as research instrumentation 
(and, under full scale conditions, modeling based on physical and 
mathematical models) -- have been formulated. That is making it 
possible to solve both many new problems in materials science 
itself and new problems in interdisciplinary fields (such as 
space biology, geophysics, and astrophysics). 


Capabilities for varying by several orders of magnitude the force 
of gravity, the most universal force in the universe, as well as 
the immense variety associated with gravity-sensitive processes 
in man’s environment and in equipment systems and with production 
and biotechnological processes, are creating new prospects for 
basic research and for practical applications. 


The Department of Problems in Machine Building, Mechanics, and 
Control Processes, together with the Institute of Problems in 
Mechanics of the Russian Academy of Sciences, held a science 
meeting called "Problems in Mechanics Under Conditions of 
Microgravity" on 11 November 1993 at the Institute of Problems in 
Mechanics of the Russian Academy of Sciences (in Moscow). 


The program for the meeting included one morning session with 
five papers given orally and an evening session that was in the 
form of a round table discussion on the topic "Models in 
Mechanics for Studying Production Processes Associated With 
Microgravity and for Posing New Problems." In the break between 
the sessions, the meeting’s participants were provided with 
additional information in the form of poster reports, 
demonstrations of computer programs, and an exhibit of the most 
recent publications in that field from academy and sector 
institutes, higher educational institutes, and aerospace design 
bureaus. 
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The papers presented theoretical and experimental work and 
methods of laboratory modeling under conditions of microgravity, 
as well as results from full-scale experiments. 


The participants from Moscow and Moscow Oblast represented the 
Russian Space Agency and the following Russian Academy of 
Sciences institutes: A. A. Blagonravov Machine Building 
Institute, A. A. Baykov Metallurgy Institute, the Solid Physics 
Institute (Chernogolovka), the M. V. Keldysh Applied Mathematics 
Institute, the Physical Chemistry Institute, the State Institute 
of Rare Metals, the Institute of Biomedical Problems, the MGU 
Institute of Mechanics, the Moscow Aviation Institute, the Moscow 
Physical-Technical Institute, the N. E. Bauman MGTU Applied 
Mathematics and Mechanics Scientific Research Institute, the . 
Moscow Aviation Institute Applied Mathematics and Electronics 
Scientific Research Institute, the Central Scientific Research 
Institute of Machine Building, General Machine Building Design 
Bureau, the NPO Energiya, the S. A. Lavochkin NPO Kompozit, the 
Science Center in Zelinograd, the Scientific Research Institute 
of Thermal Processes, and the N. Ye. Zhukovskiy Central 
Aerohydrodynamics Institute. Also taking part were staff members 
of the Physics and Energy Institute (Obninsk), the Scientific 
Research Institute for Electronics Materials (Kaluga), the NPO 
Foton (Samara), the Russian Academy of Sciences Elektropribor 
Scientific Research Institute and A. F. Ioffe Physical-Technical 
Institute (St. Petersburg), Rostov University (Rostov-na-Donu), 
the Institute of Solid Media Mechanics of the Urals Department of 
the Russian Academy of Sciences and the Perm State University 
(Perm), and the M. A. Lavrentyev Institute of Hydrodynamics of 
the Siberian Department of the Russian Academy of Sciences 
(Novosibirsk). Such a broad range of participation in the meeting 
by the various organizations points to the great amount of 
interest in our country in this area of science. 


In light of the fact that many of the problems in this area are 
based in methods associated with fluid and gas mechanics, the 
editorial board of the IZVESTIYA RAN. MECHANIKA ZHIDKOSTI I GAZA 
journal decided to publish the material presented in the science 
meeting in one of the journal’s issues for the perusal of the 
readers. Presented here are papers from the morning session, a 
number of poster reports, and the remarks made by the 
participants in the round table discussion, as well as articles 
that are related to this area of research and were submitted to 
the journal’s editor. 


At the same time, some of the work reported in this issue is not 
directly related to fluid and gas mechanics. The editorial board 
saw fit to publish them in order to give the reader a fuller idea 
of the current problems in this field of science that interfaces 
with mechanics, physics, and materials science, as well as with 
the mechanics and technology of space flight. 








Mathematical Modeling and Measurement of Microaccelerations 
Aboard the Mir Orbital Station 


957Q0020B Moscow IZVESTIYA ROSSIYSKOY AKADEMII NAUK MEKHANIKA 
GHIDKOSTI I GAZA in Russian No 5, Sep-Oct 94 [manuscript 
submitted 13 Jan 94] pp 5-14 


[Article by M. Yu. Belyayev, S. G. Zykov, S. B. Ryabukha, V. V. 
Sazonov, V. A. Sarychev, and V. M. Stazhkov, Moscow; first 
paragraph is abstract] 


[FBIS text] Causes of microaccelerations aboard artificial Earth 
satellites are examined: motion around the center of mass, 
gravity field gradient, atmospheric resistance, etc. Methods are 
described for determining the low-frequency microacceleration 
component aboard the Mir orbital station from the data of a 
three-component magnetometer and an angular velocity sensor, and 
the results of such determination are given. Also given are the 
results of microacceleration measurements made with an 
accelerometer. 


1. Microaccelerations aboard artificial Earth satellites. The 
mechanical nature of microaccelerations aboard artificial Earth 
satellites consists in the following. The orbital motion of the 
satellite is described with the equation 


R=g9(R, t) +F/m 


Here R is the geocentric radius-vector of the center of mass of 
the satellite; the dot represents the differentiation in time t; 
g(R, t) is the intensity of the Earth’s gravitational field (we 
are disregarding the gravitational influence exerted on the 
satellite by the Moon and other bodies of the solar system); m is 
the mass of the satellite; and F is the principal vector of the 
forces of nongravitational nature that are acting on the 
satellite. We fix on the satellite some point P. Its radius- 
vectors relative to the center of the Earth and the center of 
mass of the satellite we designate with, respectively, r and p: r 
= R +p. We shall label as a microacceleration at point P the 
magnitude n = g(r, t) - r, which is the difference between the 
intensity of the gravitational field at that point and the 
absolute acceleration of the latter. If at point P there is an 
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experimental body with a mass of m, << m, then the apparent force 
acting on that body is equal to m,n. Based on the correlations 
written out, the microacceleration can be represented in the form 


n=-$+g9(R+p,t) -g(R, t) -F/m 


Let the satellite be a solid. Then, designating its absolute 
angular velocity as o, we have 


b=@xp+o~x (@xp) 


With an accuracy adequate for the calculations, we may assume 
that g(R, t)=-u,R/R’®, where R* =|R| and uy, is the gravitational 
parameter of the Earth. In that context, if p = |p| << R, the 
following formula holds true 





g(R+p,t) -g(R, ¢) -2 [ : (Rp) R _ 9) 











R2 
Now the expression for n can be written in the form 
n=px@+ (exp) «a= [= ap) Eo] -< (1.1) 
Nongravitational disturbances acting on Salyut-6, Salyut-7, and Mir stations | Typical values for |F/al, 
in orbit n/s! 
Effects of control system during orbital correction 107-1 
Effects of attitude-control syste 1072-107! 
Operation of fans, compressors, and other gear 107-107 
Calisthenics of crew members 1072-107 
Aerodynamic drag 10 
Light pressure 10 
Magnetic field of Earth 10 
Micrometeorites 10 





The expression obtained explains the sources of the 
microaccelerations aboard the satellite-solid. Generallv, there 
are three sources: (1) the motion of the satellite relative to 
the center of mass — the first two terms on the right side (1.1) 
account for that motion; (2) the irregularity of the Earth’s 
gravitational field — it is described by the third term; (3) 
forces of a nongravitational nature that act on the satellite — 
the term "-F/m" describes them (such forces consist of various 
kinds of controlling factors, factors that are due to the 
movement of internal parts of the satellite, aerodynamic drag, 
light pressure, etc.). If there are simple design formulas for 
the microacceleration components stemming from the first two 
sources, deriving a clear formula for F/m requires assigning a 
number of specific mechanical models: models of the interaction 
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of the satellite and the outside medium, models of the operation 
of its controls, etc. The table gives values of |F/m| typical of 
satellites such as the Salyut-6, -7, and Mir stations for several 
disturbing influences. The following example makes it possible to 
estimate the importance of those influences. Let us examine a 
satellite in a circular orbit that is about 340 km in altitude; 
the satellite .s stationary in the orbital system of coordinates. 
Let us assume that the vector p is collinear with the vector R 
and has a length of 2 m. Then the modulus of the sum of the first 
three terms on the right side (1.1) equals 3y,9/R°® = 7-10° m/s’. 


Since it is impossible to completely eliminate microaccelerations 
aboard a satellite, analysis of the results of in-flight 
experiments requires that we have information on the magnitude of 
the microaccelerations. We can determine what the 
microaccelerations are either by using special accelerometers or, 
if we have information on the motion of the satellite, by using 
formula (1.1). Those two approaches are different from each 
other, but each has its own field of application. Accelerometer 
readings themselves do not address the sources of 
microaccelerations, and an analysis of those readings requires 
additional information. The second approach is limited to a 
satellite-solid and consists essentially in the mathematical 
modelling of microaccelerations in situations close to actual. In 
many ways, the second approach is simpler than the first. So let 
us begin with the second and demonstrate it with an example for 
determining the microaccelerations aboard the Mir orbital 
station. 


Two methods of using formula (1.1) will be examined. We shall use 
the first method only for uncontrolled motion of the station. 
With that method, the actual motion of the station around the 
center of mass is determined as a result of statistical 
processing of the data of onboard-magnetometer measurements, and 
then the microacceleration at any given point of the station is 
computed for the motion as a function of time. The solution of 
that type of problem for the Salyut-6 and -7 stations is 
described by Sarychev et al.* We shall use the second method for 
station motion that has a rather large angular velocity. With 
that method, the data of onboard measurements of angular velocity 
are smoothed with splines and are inserted into formula (1.1), in 
which the two final terms are discarded. 


2. Calculation of microaccelerations aboard the Mir orbital 
station in uncontrolled flight. This problem is of interest for 
two reasons. First, production experiments aboard the Mir station 
are better performed in uncontrolled flight, which ensures a low 
level of disturbing microaccelerations. Second, there is an 
effective methodology for determining the uncontrolled motion of 
that station from the readings of the onboard magnetometer .? 


Of the nongravitational forces acting on the station, we shall 
consider only the force associated with aerodynamic drag, which 
we approximate with the formula 
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F=mbp,vv, v=|v| 


Here b = const is the ballistic coefficient, p, is the density of 
the aerodynamic stream approaching the station, v is the velocity 
of the station relative to the Greenwich system of coordinates. 


Here we shall provide a brief description of the methodology for 
determining the motion of the station from the magnetometer 
readings and the mathematical model underlying it in the portion 
that is of interest for solving the problem under consideration. 
The station (the base module with spacecraft and other modules 
docked to it) is regarded as a solid whose center of mass 
performs a Kellerian elliptical motion around the Earth. The 
elements of that motion and the ballistic coefficient b are 
determined from trajectory measurement data. To write the 
equations of motion of the station relative to the center of 
mass, we shall introduce two right-handed Cartesian coordinate 
systems: the orbital system X,X.X, and the system formed by the 
main central inertial axes of the station x,x,x,. The axes X, and 
X, are aligned along vectors R and R x R; axis x, is close to the 
longitudinal axis of the base module and projects toward its 
equipment section; the unsymmetrical solar array of the base 
module lies in the half space x, 2 0. The matrix for the 
transition from the x,x,x, system of coordinates to the X,X,X, 


system we shall designate as 1a.jli, 5-1, @,;=COS(X;,X;) . 
The equations for the rotational motion of the station are 


written with consideration for the gravitational and restorative 
aerodynamic moments 


@, = (@,0,-Va,24,,) +X (UD; -U5P,) 


» - 1-A ’ 
aS Wri (@,@,-Va,,4,,) “i (U,D,-U,D;) 
@,=(1-A+Ap) (0,0,-Va,,a,,) tAx (U,D,-v,D,) (2.1) 


4,1 =4,20,-4,,0,-W.43,, 4,2=4,,0,-4,,0,-W,a,, 
4), =4,,0,-€,20,-W.4,,, 43, =4,,0,-4,,0,-Wa,, 


432 =4330,-€3,,0,-W.4,,, 4,3 =43,0,-2,20, -Wo a; 





Here o, and v, are the components of the vectors # and v in the 
system of coordinates x,x.,x,; I, represents the inertial moments 
of the station relative to axes x,; the parameters p, 
characterize the aerodynamic moment acting on the station; o, is 
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the modulus of the absolute angular velocity of the orbital 
system of coordinates; and E is the scaling coefficient. In the 
numerical integration of the equations in (2.1), the unit of time 
measurement is 1000 seconds; the unit of distance measurement is 
1000 km; the density of the atmosphere is calculated in kg/n,, 
according to the Model of the upper atmosphere for ballistic 
calculations. GOST 22721-77°; and E = 10’°. The missing elements 
of the a,, transition matrix are computed with the formulas 

@2, = 43;2€,, ~ €@,,;8,,, etc. 


The readings of the onboard magnetometer are interpreted in the 
so-called structural system of coordinates y,y.y,. That system is 
a right-handed Cartesian system whose axes form small angles with 
the axes of the system x,x.x,. The axis y, is parallel to the 
longitudinal axis of the base module and projects toward that 
unit’s equipment section; axi_ y, is parallel to the axis of 
rotation of the unsymmetrical solar array. The matrix for the 
transition from the system x,x.x, to the system y,y.y, we shall 


designate as Ib, li, j-1. Dj;=COS (y;,X;) . The elements of that matrix 


are expressed as a function of the angles a., 8., and y., at which 
the system y,y.y; must be turned successively around axes y,, y3;, 
and y, in order to convert it to the system x,x,.x;,. 


The values for the angles a., 8, and y. are calculated along with 
the inertial parameters A and yw from the components of the tensor 
of inertia of the station in the structural system of coordinates 
on the basis of the results of shaping that tensor into a 
diagonal. The tensor of inertia of the station is, as a rule, 
known accurately enough, and the method of assigning the 
parameters A, wu, a., B., and y. makes it possible to 
satisfactorily match the measured and the calculated components 
of the vector of the intensity of the geomagnetic field. It is, 
however, possible to refine those parameters while processing the 
measurement information along with the initial conditions of 
station motion and aerodynamic parameters p,. 


At given moments of time, measurements of the components h, (i = 
1, 2, 3) of the vector H of the intensity of the geomagnetic 
field in the structural system of coordinates are taken from the 
magnetometer. Usually, not all accessible measurements are 
included in the processing — roughly only every fourth to eighth 
measurement is included, which means that the interval of time 
between neighboring measurements of the same component is several 
minutes. Such sampling is adequate for determining the rotational 
motion of the station with an angular velocity of about 10° c”. 
The aggregate of measurements selected for processing we shall 
designate as 


ti". hj”, (k=1,2,...,N,; i=1,2,3) (2.2) 











Here h:*’ is the result of the measurement of the component h, of 


+ 


(k) 


the vector H at the moment in time ¢;”’, t/*<t/***). The beginning 


and the end of the time interval being processed we define with 


: ; (Ni) 
the correlations t,=mint;*’, t,=max ¢t; (i=1,2,3). As a rule, 


t, - t, is less than or approximately equal to 120 min, and N, = 
30-50. 


When the motion of the center of mass of the station are known, 
the components of the vector H in the orbital system of 
coordinates can be calculated from well-known formulas. 
Determining the rotational motion of the station in the segment 
t. < t < t, consists in seeking out the parameters of the 
mathematical model and the solution of system (2.1) that best 
allow, in the context of the least squares method, matching the 
measured and the calculated values of the components of H. 
Processing of the measurements is broken down into several 
stages. 


First, on segment t, < t < t,, cubic splines H,(t) (i = 1,2,3) are 
constructed that express the design time dependence of the 
components of vector H in the orbital system of coordinates. A 
grid with a spacing of about 2 minutes is used, and the values 
for the geomagnetic field at the nodes of the grid are computed 
with the Golovkov and Kolomiytseva model.* Then the statistical 
processing itself of the measurement information of (2.2) is 
performed. The functional of the least squares method is formed 
on the basis of the following hypothesis: the systematic errors 
in the measurements of the h, component are identical and are 
equal to A, (i = 1,2,3); random errors in the measurements of all 
components are independent and have an identical normal 
distribution with a zero mean value and a standard deviation o. 
That functional take the form 


3. ON; 3 
®=> > (hj -h? (7%) 17, h?(t) =Ay+ YO bj ,;a,. (t) A, (t) 


171 k=1 j, 11 


The problem consists in minimizing that functional on the basis 
of the initial conditions of solution of system (2.1) that are 
assigned at point t, and on the basis of the parameters p, and A. 
(and, possibly, A, wu, a., B., and y¥.). 


After @ is minimized, the solution that is found for system (2.1) 
can be used to calculate the microaccelerations on segment t, < t 
< Cee 


As an example of that, let us cite the results of the calculation 
of microaccelerations for a measurement interval for 16/17 March 
1993. In that interval, t, = 23 hr 37 min 12 sec and t, = 1 hr 45 
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min 5 sec, official Moscow time; NW, = 50, N, = 53, N, = 52; after 
refining all the above parameters for o, an estimate of o = 7307 
was obtained. The precision of the match of the measured and the 
calculated values for h, is illustrated in Fig. 1. In that 


figure, the solid lines depict plots of the functions A°(t) that 


correspond to the point at which @ is at its minimum, and the 
circles indicate the measurements included in the processing. 
Figure 2 presents the plots of the time dependence of the 
components n, (i = 1,2,3) of microacceleration n in the system 
y.:y2y; that are conctructed for the vector p, which, in that 
system, has the comporents (-9 m, 12 m, 1m). In this case, the 
vector p was chosen arbitrarily. 
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A microacceleration calculated with that method is a background 
function of time. In reality, fluctuations and individual peaks 
resulting from jolts, station hull vibrations, and other sources 
are superimposed on it. If the station is regarded as é solid - 
and that is precisely the assumption under which formula (1.1) 
was derived — it is impossible to allow for such effects. For 
that reason, the method provides a sufficiencly accurate idea of 
only the low-frequency microacceleration component (frequencies 
less than or approximately equal to 10° Hz). 


3. Determination of microaccelerations from the data of onboard 
measurements of angular velocity. Aboard the station, there are 
angular-velocity sensors that make it possible to measure, at 
given moments in time, the components Q, (i = 1,2,3) of the 
vector # in the structural system of coordinates. Thos<« 
measurements can be made both during uncontrolled flight and when 
the station is being acted on by controlling forces and moments. 
The aggregate of the measurements of e that are performed in sone 
time interval we shall designate as 


c/* ,Q:” (k=1,2,...,N,; i=1,2,3) (3.1) 


Here Q* represents the results of the measurement of the 


components f, at the moment in time t/*’, ¢/* <el*? 


_ - 


As a 


rule, 7 -ti* = 20-60 s and NW, = 10-50. The solution of the 


following problem is ued in the processing of the data of the 
measurements of (3.1). 


For moments in time ¢, (k=1,...,N), t,<t,.,, let us assume the 
approximate values of x,=f(t,)+e, of some smooth function f(t) and 


rms estimates of $x, of errors e,. That function needs to be 
reconstructed for the segment t, < t < t,. In a paper by Reinsch’, 


the search for f(t) under the assumption that that function is 


differentiated twice without interruption boils down to a 
solution of a covariational problem 


r 
tJ 


N 

: : X,~f£ (,) 

= 2/ ~mi ( k k 
J-f? t) dt min, ) 


5x, 





)*<S (3.2) 


Here S is an assigned positive number. The solution of the 
problem (3.2) is a cubic spline. Reinsch’® presents a program that 
comp tes the coefficient of that spline for S, t,, x,, 5x, (k = 
1,...,N). That program is used for plotting the 
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splines Q,(t) that smooth the measurements of (3.1), i.e., that 


satisfy the conditions Q,(t{") is approximately equal to Q}*’. 


In the context of that covariational problem, the parameter S is 
undetermined. There are, however, a number of statistical 
formulas for choosing suitable values for it based on the 
magnitudes of t,, X,, and 6x,. In smoothing the measurements of 
(3.1), the following technique was used (compare Wahba‘). Let us 
assume that s 2 0, j is natural number, and f, (t, s) furnishes 
the minimum to the functional J in (3.2) under the condition that 


N 


» ( th) yags 
“1 Xx 


ke) 





That function is a cubic spline that can be plotted with 
Reinsch’s program.* Let us take A to be greater than 0 and p to 
be a natural number from the interval (1, N/2). We shall assume 
that 





Np*l =x. -f,(t., 8) 
= | Jj & 2 
® (Ss) » ( 3 ; ) 


mooguis @(10-°A) , b=0,1,... 


It is assumed that S = 107AN. 


Examples of the splines that smooth the measurements of (3.1) are 
given as solid lines in Fig. 3. The circles in that figure 
indicate measurement data. The measurements were made during a 
Eulerian turn of the station about an axis that was stationary in 
the structural system of coordinates and in absolute space. Here 
N, = N, = 17, Nz = 13, and the moment t = 0 corresponds to 19 hr 1 
min 1 sec, official Moscow time on 4 December 1993. The splines 
were plotted for 6x, =1, p = 3, A= 1, and B = 3 for all three Q, 
components. An example of the use of those splines for 
microacceleration calculation is given in Fig. 4. The 
microacceleration is calculated with the 


formula n=px@+(@xp) x@ for the vector p = (-9 m, 12 m, 1 m). 


That formula is applicable if the sum of the first two terms on 
the right side of (1.1) exceed in absolute value the sum of the 
two final terms. In that case, the condition is fulfilled with 
stress. 


This method of microacceleration calculation — just like the 
preceding one — makes it possible to actually find only the low- 
frequency component. The upper boundary of the computed function 
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of n(t), however, is greater than that of the first method by at 


least an order of magnitude. 
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4. Measurement of microaccelerations with an accelerometer. At 
present, measurements of microaccelerations aboard the Mir 
orbital station are performed with two sets of Mikroakselerometr 
gear developed by CNES (France). The gear consists of an 
electronic unit with a video recorder, a three-component unit of 
sensors with three types of mechanical interfaces, and a video 
camera that enables monitoring of the location where the sensor 
unit is installed and of the orientation of the axes of their 
sensitivity in the structural system of coordinates. The 
Mikroakselerometr enables measurement and recording of 
microaccelerations with an error of around 1 percent in a range 
of 10°°-1 m/s’ and a frequency band of 0.1-400 Hz. Measurement 
results are recorded on video cassettes that are returned to 
Earth for detailed analysis. In addition, it is possible to relay 
the measurement information via video channel for day-to-day 
monitoring of the level of microaccelerations. 


Statistical and spectral analysis of the measurement data show 
that the level of microaccelerations on the Mir station is 
substantially higher than the levels aboard Salyut-6 and -7. The 
average values of |n,| depend on the time and place of the 
measurements and lie with a range of 3-10°7~15-10°? m/s?; the 
average values of |n,| and |n,| can reach 1 m/s?. Most of the 
strength of the microaccelerations is concentrated in a frequency 
range of 20-50 Hz. In the station’s base module, the principal 
source of disturbances is the BKV-3 compressor. The angular 
velocity of its proper rotation is 24.4 Hz. That frequency and 
multiples of it are quite clear in the spectrum of the function 
|n(t)|. In the vicinity of the central post of the base module, 
effective microacceleration-component amplitudes of 0.66 m/s? and 
0.21 m/s’? correspond to spectral peaks with frequencies of 48.8 
Hz and 146.4 Hz. When the compressor is shut off, the level of 
microaccelerations in the base module drops by an order of 
magnitude. 


The gyrodynes are another substantial source of disturbances. 
When the gyrodynes are in normal operating modes, the 
Mikroakselerometr registers microacceleration levels with 
effective amplitudes of about 10°? m/s?. The spectrum of 
microaccelerations shows peaks at frequencies of 150-160 Hz and 
300-320 Hz that are due to the proper rotation of the rotors of 
the gyrodynes with angular velocities of around 10,000 rpm. The 
amplitudes of the spectral components corresponding to those 
peaks are as high as 0.02 m/s’. Along with the high-frequency 
components in the microaccelerations, there are also low- 
frequency components with frequencies of 25 Hz and 45 Hz and 
amplitudes of 0.01-0.02 m/s’. But during dynamic operations 
(turns, unloading, etc.), the amplitudes of the low-frequency 
components of the microaccelerations rise abruptly, roughly by an 
order of magnitude. In addition, those amplitudes also rise 
abruptly when certain other station systems whose natural 
frequencies are close to the typical frequencies of the gyrodynes 
are switched on. 
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The Mikroakselerometr gear has enabled measurement of the level 
of microaccelerations due to elastic vibrations of the station 
hull during calisthenics performed by crew members. Such 
measurement was the purpose, on 22 December 1992, of the Resonans 
experiment, in which elastic vibrations of the hull were produced 
by a cosmonaut running on a treadmill. During the experiment, 
many of the station’s systems that make a substantial 
contribution to the level of microaccelerations aboard the 
station were switched off. The Mikroakselerometr was positioned 
next to the Kseniya telescope in the Kristall module. Before the 
vibrations were produced, the level of microaccelerations |n| 
=0.007 m/s’ was recorded; and 16 percent, 20 percent, and 10 
percent, respectively, of the power of the function |n(t) | were 
concentrated in the frequency ranges of 0.10-3.61 Hz, 8.59-13.48 
Hz, and 50.78—-55.66 Hz. After the vibrations were excited, the 
level of microaccelerations rose to 0.026 m/s’, and 18 percent 
and 77 percent of the power of the function |n(t)| were 
concentrated in the frequency ranges of 0.10-2.34 Hz and 
9.67-13.96 Hz. The spectra of that function before and after 
excitation of the vibrations are given on the upper and lower 
graphs, respectively, of Fig. 5. The graphs are plotted on a 
logarithmic scale; frequency in hertz is plotted along the x- 
axis, amplitude in m/s along the y-axis. The upper graph is a 
plot of measurements performed for a time interval of 60 seconds; 
the lower graph, of measurements for a time interval of 107 
seconds. In both cases, the segments of time between the two 
consecutive measurements was 0.005 s. It follows from the data 
that the calisthenics and work performed by the crew members 
makes a contribution to the level of microaccelerations aboard 
the station that is roughly an order of magnitude smaller than 
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the contribution made by the various control and life-support 
systems. 


The work in this paper was performed with the support of the 
Russian Fund for Basic Research (project code 93-013-16249). 
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System for Protecting Onboard Production Equipment from 
Vibrations 


957Q0020C Moscow IZVESTIYA ROSSIYSKOY AKADEMII NAUK MEKHANIKA 
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[Article by F. K. Kalashnikov, V. L. Levtov, L. V. Leskov, V. V. 
Romanov; first paragraph is abstract; UDC 531.3:534.1:629.7] 


[FBIS text] Various methods of protecting onboard production 
equipment are analyzed in this paper, as are the results of the 
Vibratsiya experiment, which was conducted aboard the Mir manned 
orbital station. 


1. Effect of microaccelerations on production processes. Analysis 
of the results of production experiments conducted aboard manned 
orbital stations and unmanned space vehicles has shown that one 
of the chief factors preventing the potential advantages of the 
production of materials in weightlessness from being fully 
realized consists in quasiconstant and variable 
microaccelerations.’ For the in-orbit production of semiconductor 
materials, those advantages consist of the following: higher 
degree of structural perfection of the single crystals, 
composition close to stoichiometric, higher degree of purity, 
higher degree of uniformity of electrophysical parameters. For 
the production of ultrapure biological preparations via 
electrophoresis, the advantages are a higher resolution and, 
accordingly, higher degree of refinement, as well as higher 
productivity for the process itself. Similar potential advaritages 
exist for other types of materials. 


The actual dynamics aboard an unmanned vehicle and particularly 
aboard orbital stations make it possible to realize those 
advantages to only a small extent. When samples of silicon and 
gallium arsenide are grown aboard the Mir orbital station, for 
example, the melt separated from the growing crystal in some 
instances. According to estimates, the magnitude of 
microaccelerations was 6-107? g, (g, = 9.8 m/s’). If there is to be 
high-quality growth of single crystals with a diameter of 20-50 
mm, the magnitude of microaccelerations cannot exceed 10° g,. 
That pertains both to the quasiconstant component and the 
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vibrational accelerations in the most hazardous frequency range 
of 0-100 Hz. 


Estimates of the quasiconstant component of microaccelerations 
due to the operation of control systems yield a magnitude of 10° 
2-10" g, for the Mir station; that figure is as high as 0.1 g, for 
maneuvers perfcrmed for correcting station trajectory and 
attitude. Measurements of the vibrational background in the 
frequency range of 20-50 Hz yield a magnitude of 
microaccelerations of 3-10°°=—..5+-107? g,; in individual modes, the 
figure is as high as 0.1 g,. That means that, if special measures 
are not taken to protect the production equipment from vibration, 
only certain production experiments of a basic nature should be 
performed aboard the orbital station, and full-scale processing 
of engineering production processes - chief among them the 
production of semiconductor materials — are not possible. 


As for problems of space-based biotechnology, similar 
quantitative estimates of the allowable level of 
microaccelerations have not yet been made. We know, however, that 
when processes such as electrophoresis are used, the generation 
of the Joule effect results in a substantial elevation of the 
temperature of the biological preparations.” Since the Rayleigh 
and Grashof numbers, which characterize the intensity of gravity- 
induced convective flows in a fluid, are proportional to the 
third degree of the linear dimension of the zone of flow, there 
is reason to expect that, in that case, the microaccelerations 
will lead to negative effects. 


What has been said above means that one of the central problems 
on the road to realizing the potential advantages of space-based 
materials production consists in developing special systems for 
protecting production equipment carried aboard space vehicles 
from vibrations. 


2. Methods of protecting against vibrations. After the 
performance of the first experiments that made it possible to 
demonstrate the effectiveness of protecting equipment from 
dynamic loads, the need arose to systematically study various 
versions of the solution of the problem. Let us consider which 
methods can be used for that and what their advantages and 
disadvantages are. 


We are aware of various kinds of devices (vibration insulators) 
mounted between the source of the vibration and the object being 
insulated. Various elastic elements can serve as a vibration 
insulator: a spring, a cable, rubber, etc. If the response of the 
elastic elements is taken to be static, then we get a static 
force curve of F = cx, with the insulators divided into linear or 
nonlinear, depending on the type of curve. 


Note that when an elastic element is mounted with minimal 
rigidity, the dimensions of the vibration insulator, which are 
limited by the operating conditions of the object being 
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protected, grow. The noted drawbacks, along with others that are 
inherent to linear elastic elements, force us to go to the use of 
nonlinear systems of vibration protection. 


This method of vibration protection is a passive system. Active 
vibration-protection systems are created to protect against low- 
frequency vibrations and variable vibrational spectra. 
Characteristic of an active vibration-protection system is the 
presence of an additional source of energy and controlling 
devices. The chief advantage offered by active vibration- 
protection systems is a hign supporting capability with low 
dynamic rigidity. That, however, is achieved by making the design 
of the vibration-protection system considerably more complicated. 
But a given supporting capability, with low rigidity, can be 
effected with a simpler method — by combining elastic elements 
that are linked in a specific way. That class of vibration- 
protection systems -— called vibration-protection systems with 
quasizero rigidity — will be described below. 


3. Passive vibration-protection systems with quasizero rigidity. 
One way of getting a segment of quasizero rigidity on the 
performance curve of a vibration-protection system is through a 
parallel hookup of a rigidity corrector — i.e., a device that has 
a "negative rigidity" segment in its force curve - to a linear 
vibration-protection element. That could be, for example, a 
sinusoidally bent rod, precompressed from the ends to loss of 
stability; when transverse force is applied to it, it "cracks" in 
the direction of the action of the force, i.e., it has "negative" 
rigidity. 


A fundamental property of systems with quasizero rigidity is that 
they, being tuned to a zero supporting capability (and that’s the 
main condition of weightlessness), have the same force curve that 
they have on the ground. The only difference is that now the 
performance segment of the quasizero rigidity on the force curve 
either intersects the x-axis at a small angle or lies on the x- 
axis. That indicates the possibility of high-quality vibration 
protection in one or several directions in space. Thus, if a 
rigidity corrector is attached to the elastic element, such a 
system in weightlessness affords high-quality protection of the 
object against vibration. The shape of the linear element 
(cylindrical, disk, spiral, etc.) is of no importance whatsoever; 
the corrector, by dint of its own properties, adjusts the 
rigidity all the way to zero and affords low-frequency protection 
in weightlessness. Even if the supporting spring is nonlinear, 
there will be correction (albeit of a lower quality). 


In a vibration-protection suspension designed in the form of a 
"spring cross," there are no correctors as individual elements, 
and the force curve with a segment of quasizero rigidity is 
achieved by selecting the proper correlation of linear elastic 
elements (usually standard cylindrical springs). The advantage of 
that method of suspension is simplicity and producibility. 
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Fig. 1 Fig. 2 


Figure 1 presents one of three shock absorbers with such a 
suspension. It is made in the shape of two yokes, 1 and 3, which 
are set at a 90° angle to each other. Each yoke is closed by a 
pair of series-connected springs 2. All four springs of the shock 
absorber are soldered together at a single point A. Yoke 1 is 
secured to the object being protected, and yoke 3 is secured to 
the base. Mounting the object on those shock absorbers ensures 
high-quality protection of the object against vibration. 


The force curve of the system in shown in Figure 2 for a rigidity 
of vertical springs of 10 g/cm and for a maximum deviation of the 
system from the position of equilibrium of X, = 10 mm. For 
purposes of comparison, the linear curve of a standard spring 
with a rigidity of 0.25 g/cm is shown with the broken line, i.e., 
if a spring 40 times more rigid were used in that arrangement, 
the rigidity of the suspension would be less. That produces a 
positive effect on the attenuation of higher resonances of 
springs by means of raising their frequency. It can be seen from 
the graph that the suspension has a nonlinearity and that there 
is a segment of quasizero rigidity in its force curve. Such a 
suspension can be used for small masses (about 1 kg) and large 
masses. 


That suspension does not require tuning or adjustment of the 
corrector, which is a weak point in other arrangements of the 
same class. It is simple and reproducible, requires no current 
supply or control system, is explosion-safe, does not generate 
toxic substances, and creates no electrical or magnetic fields. 
The suspension makes it possible to easily solve problems of 
current supply to the unit via the use of the shock absorber 
springs — which are insulated and coated with copper ahead of 
time — as a conductor. 


It is possible to use the cables currently existing on the unit 
(after some modification, perhaps) as the main elastic or damping 
elements of the design without any substantial reduction in the 
quality of vibration protection and, possibly, with an 
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improvement in such protection. 
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Fig. 3 Fig. 4 


Figure 3 presents a different type of vibration-protection 
suspension with elastic elements in the form of cable arcs 2, 
which are used to secure the object being protected to the hull. 
For lowering the natural frequency, it is equipped with a spatial 
corrector 3, which takes the form of plates or rods precompressed 
from the ends to loss of stability. The corrector makes it 
possible to reduce the rigidity of the structure. Overall, that 
design has a number of merits: low rigidity, excellent damping of 
the main and the higher resonances, applicability for small 
masses (about 1 kg) and large masses. 


Another example of a nonlinear cable suspension is a suspension 
with a needle rigidity corrector. The supporting elastic element 
takes the form of a cantilever made of steel cord or cable i 
(Fig. 4). The suspension has a linear force curve and, 
consequently, has all the advantages and disadvantages of linear 
vibration-protection systems. To lower the natural frequency, the 
suspension is equipped with a rigidity corrector in which the 
needle 3 is thrust between a ruby prism 2, which is rigidly fixed 
to the free end of the cantilever 1, and a leaf spring 4. By 
tightening the needle with the leaf spring, one can achieve an 
unstable situation with the needle-—cable system, and the system 
acquires the property of "cracking," i.e., a negative rigidity. 
Sy adjusting (5) the tightness of the leaf spring, one can 
achieve the desired force curve with a segment of quasizero 
rigidity. The restoring force of the suspension shown in Fig. 4 
lies in the I-I plane. By connecting two modules in the form of 
two yokes set at 90° to each other (as in the "spring cross"), we 
obtain a spatial shock absorber with quasizero rigidity, and its 
oo. 1% of operation is similar to that of the "spring cross" 
(Fig. 5). 


That suspension has a number advantages. First, the use of a 
cable spun from fine steel wires as the elastic element made it 
possible, as a result of the internal friction in those wires, to 
virtually eliminate the higher resonances (that is also 
facilitated by the simple shape of the elastic element - a 
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cantilevered beam), as well as attenuate the coefficient of 
amplification on the main resonance (for cable suspensions, it is 
2-3). Second, the elastic element or the cable cores can be used 
as conductors if they are copper-plated and insulated. And third, 
the use of a needle corrector whose pointed end rests in the ruby 
prism makes it possible to lower the force of dry friction and 
raise the sensitivity of the system, i.e., make the suspension 
efficient in protecting equipment of small masses of around 1 kg, 
when high sensitivity is particularly necessary. 














Fig. 5 Fig. 6 


For protecting units with a mass of around 1,000 kg from 
vibration, when a fluid feed is requires, one can utilize a 
Similar design, using the pipe itself 1, bent in the shape of the 
letter [T as the elastic element (Fig. 6); one end of the pipe is 
secured to a flange of the base 4, the other end, to a flange on 
the equipment 5. At the bend in the pipe, two needle correctors 2 
are mounted, making it possible to lower the rigidity of the 
supporting pipes. The correctors are tightened by a vertical leaf 
spring 3 and a horizontal spring 6, which is filled with damping 
material to lower the natural resonances. The higher resonances 
of the elastic element itself are damped by the fluid in the 
pipe, whereas the main resonance of the elastic element can be 
damped by placing electric cables on the pipe that are secured 
with clamps. 


4. Magnetic vibration-protection systems. The action of systems 
of this type of vibration protection is based on the property of 
magnetic fluids to be attracted to a magnet. The force of the 
magnet’s attraction per unit volume of magnetic fluid is equal to 
the product of the fluid’s magnetization and the magnetic field 
gradient and is aligned along the axis of that gradient. 
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A force of 0.5 N acts on 1 cm of standard magnetite-based 
magnetic fluid with C = 0.2 (C is the concentration of magnetite 
in colloid) in a 3-kE field with grad H = 1 kE/cm; that force is 
50 times greater than the weight of 2 cm’® fluid. For that reason, 
nonmagnetic bodies easily come to the surface in a magnetic fluid 
in a longitudinal magnetic field, with the gradient aligned along 
active acceleration. Such a field can be created with either 
permanent magnets or an electromagnetic solenoid. 


If a space vehicle is carrying an instrument that has a mass of 
100 kg and needs to be protected against vibrations with an 
amplitude on the order of 5-10% g,, then the amplification 
developed by the damping system must be around 0.5 N. When g = 
10°? g, and f = 1-1000 Hz, the depth of displacement is 0.5-10 
cn. 


It follows that if in the high-frequency region, magnetic 
vibration-protection systems can be extremely effective, then the 
use of such systems for frequencies on the order of 1 Hz or less 
runs into fairly large difficulties. Rather than permanent 
magnets, solenoid coils are probably needed through whose 
windings a 300-1000 A current will have to be passed to ensure 
the magnetic characteristics alluded to above. 


5. Choice of the base version of suspension. Summing up a 
comparative analysis of various vibration-protection systems, we 
shall formulate the basic requirements that must be satisfied by 
the base version of the suspension intended for testing aboard 
the Mir orbital station: (1) natural frequency of no more than 
0.2-—0.3 Hz; (2) minimal dimensions and weight; (3) simplicity and 
reliability of design; (4) effective in three station axes; (5) 
absence of autonomous control system; (6) capability of 
electrical feed and measurement wiring for equipment without 
appreciable worsening of vibration protection; (7) low energy 
consumption or absence of energy consumption; (8) fire- and 
explosion-resistant, and absence of toxic releases and 
appreciable electrical or magnetic noise. 


The requirements of simplicity and reliability are best met by 
systems with linear elastic elements. But when an attempt is made 
to use them to meet the requirement of low natural frequency, the 
dimensions and, as a result, the weight of the suspension grow 
substantially. For that reason, we must abandon the use of linear 
vibration-protection systems. 


Active vibration-protection systems are very effective, have 
virtually zero rigidity, and are totally free of resonances. They 
require a great deal of current, however, as well as precision 
radio engineering equipment. They are also difficult to tune, 
manufacture, and operate. For those reasons, they are not 
suitable for the orbital station currently being used. 


Nonlinear vibration-protection systems with quasizero rigidity 
best satisfy the aggregate of requirements for the base version 
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of the suspension. They are slightly worse than linear 
suspensions in terms of simplicity, but are considerably better 
in terms of satisfying the requirement of low natural frequency, 
with insignificant dimensions and weight. Moreover, they are 
resistant to the main resonance and, because of the use of 
considerably more rigid main elastic elements, have much greater 
mechanical durability. Such systems do not require current 
consumption or controlling gear. 


Thus, a nonlinear vibration-protection system with quasizero 

rigidity is suggested for low-frequency spatial protecting the 

production equipment operating aboard the Mir orbital station " 
from vibration. 


Of the designs presented in this paper for nonlinear vibration- 
protection suspensions with quasizero rigidity, the design that 
most fully satisfies the aggregate of requirements for the base 
version is the "spring cross" vibration-protection system 
described above. 


It is also pertinent to note that the known advances abroad 
either are at the stage of theoretical analysis or have been 
built as laboratory mockups; but not a single one has been 
brought to the stage of full-scale testing. In addition, judging 
from the publications available, those advances will be very 
expensive because the principles underlying their design require 
the use of highly precise technology and microprocessors. 


6. Results of the Vibratsiya experiment. One of the experiments 
was performed on 6 March 1991 (orbit No. 910) aboard the Mir 
orbital station for the purpose of measuring the principal 
characteristics of the vibration-protection suspension of the 
Vibratsiya module, which was a mass-and-dimensional model of the 
suspension planned for use to secure production units aboard 
various classes of space vehicles and was a "spring cross" (the 
mass of the object being protected was 1.5 kg). 


The Vibratsiya module was secured to the ceiling over the central 
station of the complex by elastic belts under tension. Sensors of 
the IV-1 vibrational acceleration meter were mounted on the 
vibration-protected platform of the module. An ope~ator using 
mechanical clamps either secured the platform to the hull ina 
rigid mechanical connection or "incorporated" vibration 
insulators, and the platform absorbed the station vibrational 
background reduced by the insulators. 


The degree of attenuation of vibrational accelerations on the 
platform had to be registered at various frequencies, as did the 
time it took the platform to return to the position of 
equilibrium after a short-term disturbance induced manually by an 
operator. The table presents the results of the measurement of 
the vibrational accelerations for a platform rigidly secured to 
the base and for a situation in which there is a vibration- 
protection suspension between the object being protected and the 
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base. The scale of the vibrational accelerations cited in the 
table is 10° g, (where g, = 9.81 m/s’? is the free-fall 
acceleration near the Earth’s surface). 


As can be seen from the table, before the incorporation of the 
vibration insulation, the vibration background on the object 
being protected was 2-10°-5-10% g,, with peaks of up to 2.6-10~ 
g. in the frequency range of 40-60 Hz. 
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The results of the measurement of the vibrational accelerations 
after the use of the vibration insulators show that the level of 
vibrations was reduced substantially along all the axes for 

frequencies above 1 Hz, and that included "cutting" all the peaks 
that occurred in the initial spectrum and lowering the background 
level by 20- to 50-fold, i.e., to 10* g,. 


The vibration -nsulators of the Vibratsiya module have a natural 
frequency of 0.5 Hz and, consequently, should ensure the 
reduction of the level of vibrational accelerations by an order 
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of magnitude, beginning at 1.5 Hz, which is confirmed by the 
experiment. 


The advances of the authors make it possible to create vibration 
insulators with natural frequencies of 0.05-0.08 Hz and, 
accordingly, a lower boundary for effective operation at a level 
of 0.2-0.3 Hz; theoretically, the upper boundary for effective 
operation does not exist for such insulators. Of course, in 
practice, instances are encountered in which nonmain (additional ) 
resonances due to the natural frequencies of the structural 
elements of the suspension appear in the beyond-resonance region. 
But the expected frequencies of the additional resonances are 
sufficiently great, and controlling those phenomena does not 
present any particular problems. 
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Round Table Discussion 
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[Article consists of round table discussion moderated by 
Academician G. G. Chernyy] 


[FBIS text] The round table discussion on the topic "Models of 
Mechanics for Studying Production Processes That Are Performed in 
Microgravity, and New Problems" was moderated by Academician G. 
G. Chernyy. 


Below is a brief description of the comments made in that 
discussion. 


In his comments, Academician V. S. Avduyevskiy (deputy chairman 
of the Russian Space Agency’s Problem Council for the Production 
of Materials With Improved Properties in Weightlessness) noted 
that systematic studies have been under way in our country since 
the mid-1970s in the field of hydromechanics and heat-and-mass 
transfer in weightlessness, and the results of those studies were 
discussed at all-Union seminars. The first in the series of those 
seminars was held at the Russian Academy of Sciences Institute of 
the Problems of Mechanics in December 1979; the second, in 1981, 
in Perm; the third, together with Yu. A. Osipyan and his solid- 
state physics specialists in Chernogolovka; the fourth, in 1987, 
in Novosibirsk, together with the Institute of Continuum 
Mechanics of the Ural Department of the USSR Academy of Sciences. 
The fifth, and last, all-Union seminar — held aboard a liner 
traveling from Perm to Moscow — was also known as the First 
International Symposium on Hydromechanics and Heat-and-Mass 
Transfer and gathered a broad range of specialists from the 
European Space Agency. The proceedings of all those seminars have 
been published, and the last one came out in English, in a large 
volume of papers written jointly with foreign authors. It is 
possible that the next in the series of those symposiums, which 
are now international, will take place in 1995. 


At present, the activities of our microgravity specialists are 
being coordinated by the Russian Space Agency’s Problem Council, 
as well as by the microgravity section that was created in the 
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Russian Academy of Sciences Space Council. Thus, that area, like 
the area of mechanics, has already been set up with its 
applications. 


A great deal of work is being done in the field of the 
hydromechanics of weightlessness and its applications. Certain 
areas that came about when that field of science was just 
beginning to develop intensely, i.e., about 20 years ago, have 
already lost their ties with the immediate applications and are 
serving the needs of science itself, which is also alright. At 
regularly held international conferences, one can see a multitude 
of papers on all the sections of hydromechanics that belong to 
that area (cap’llarity, convection, particle motion). 


But despite the very broad front of work being done in our 
country and abroad, there is no direct answer to the question, 
How does the change in the behavior of a fluid in weightlessness 
affect the quality of crystals being grown, and how does one find 
the optimal conditions for crystal growth, taking into 
consideration the multitude of other active factors? 


As it turns out, for example, vibrations aboard an orbital 
station play a substantial role, as a result of which crystals 
grown in space can be even worse than those grown on the ground. 
At the same time, work is being done in which a magnetic field is 
being used to suppress instabilities of convective flow and 
turbulence that take place on the ground. But here, there are 
also side effects, which were a focus of consideration at our 
seminars. 


Thus, in weightlessness, there are still no profound scientific 
results that could identify the features of phase 
transition/crystallization at the solid/liquid interface (melt). 
In recent years, however, a good many empirical facts have been 
collected, and the work being done in that area should be forces 
and combined into a single science program whose purpose could be 
to systematize those facts and identify the active factors in 
order to evaluate how each of them works. It needs to be realized 
that many of those factors are studied in isolation from actual 
crystal growth (and that includes important factors like the rate 
of crystallization), and diagnostic systems need to be developed 
in order to assess how, in actuality, a crystal grows. 
Furthermore, that question needs to be taken up by solid-state 
physics specialists (who are working in a section of the Russian 
Academy of Sciences Space Science Council under Yu. A. Osipyan) 
in order to determine what the consequences are of all those 
processes in a solid crystal and what the most important 
characteristics of such a crystal are. Thus, success in the 
solution of those problems will come in the consolidation of 
efforts of those groups of specialists. 


The Russian Space Agency’s Science Council is coordinating that 
work by setting up a general program of research. The council’s 
sphere of activity also includes aspects of international 
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collaboration, which is particularly beneficial in that field. We 
have good contacts with the European Space Agency, whose 
representative H. Walter took an active part in setting up the 
above-mentioned international symposium in July 1991. The 
European Space Agency has plans for continued collaboration in 
the field of microgravity in the context of manned flights to the 
Mir station. The prospects of collaboration with the United 
States are quite broad in connection with the agreement on the 
Mir/Space Shuttle dockings (there will be about 10). 


In conclusion, V. S. Avduyevskiy noted that, on a considerably 
smaller scale, we are doing work in biotechnology: we are 
producing new medicinal preparations via the use of a technology 
for growing protein crystals and separating substances in 
weightlessness by means of, for example, electrophoresis. There 
may be very immense possibilities here, but for now they are 
being utilized on a small scale. Unlike in space-based materials 
science, very little has been done in this field in terms of 
producing models of mechanics, and the work is of a very 
empirical nature. 


Academician Yu. A. Osipyan (chairman of the Space-based Materials 
Science Section of the Space Science Council of the Russian 
Academy of Sciences) noted in his comments that, at present, the 
science activity in the field of microgravity in our country is 
being coordinated by two science councils, one of which is set up 
in the Russian Space Agency, the other of which is a section of 
the Russian Academy of Sciences Space Council. 


In assessing the experience we have garnered and the results of 
the work that has been done since, roughly, 1976, Yu. A. Osipyan 
focused on the central problem —- the possibilities associated 
with producing better materials in weightlessness. 


Despite isolated cases of positive results, the general answer to 
the question of whether we can grow crystals of higher quality in 
weightlessness is no, because sometimes the crystals grown are of 
even worse quality. The only soluticn to the problem is to 
perform further systematic studies of the production of materials 
in weightlessness. 


One of the crucial factors active in weightlessness is the 
attenuation of convection that is due to mass forces. But 
replacing it is surface convection (Marangony convection). Thus, 
absolutely new factors are appearing that are having a 
substantial influence on the process associated with the 
formation of solids. Here solid-state physics, which is studying 
that state carefully, is not working. Research is needed in fluid 
mechanics, as is research involving the processes that occur at 
the crystallization front. Thus, new models need to be built and 
studied carefully. 


The fluid and gas mechanics specialists have the final word. That 
is why the papers that were presented at this session are 
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important for the response to the questions posed. That response 
can be given in the successive study of all the factors mentioned 
in isolation. 


Our country has performed many so-called production experiments 
in space. The Ministry of General Machine Building had most of 
the influence in the decision making regarding the conduct of 
those experiments, in the lead of which was the orbital station. 
That situation, in which crucial decisions were made without the 
input of professionals and the results could not be criticized, 
existed for a rather long time. 


When the issue of setting up the Russian Space Agency came up, it 
was decided that we need to have, as do the United States and 
Europe, commissions that would be independent of the Space Agency 
and not interested in ownership, but rather would provide 
consultation. The Russian Space Agency Science Council, which has 
now been created, plans the overall program of operations. 


For a long time, there has been lead Academy of Sciences 
institution in that field. The Space Research Institute 
laboratory headed by L. L. Regel used to play what was primarily 
only a dispatcher role. The actual research was done by the 
institutes: the Institute of Crystallography, the A. A. Baykov 
Institute of Metallurgy, the A. F. Ioffe Physical-Technical 
Institute, the State Institute of Rare Metals, the State Optics 
Institute, and a number of others. But now the laboratory has 
been reorganized. We need a science center that both would do 
science and would coordinate research. The director of the Kaluga 
branch of the Institute of Crystallization, Prof. B. G. Zakharov, 
has undertaken that work, and we need to help hin. 


As for the problem area itself, the question of the importance of 
research or production in microgravity is not clear. At present, 
the predominate trend is toward the production of semiconductor 
materials, a trend that got its start about 20 years ago with 
experiments in crystallization of germanium with traces of 
alloys. Experiments are being performed in the growth of single 
crystals of gallium arsenide, but there are no grounds to presume 
an advantage to growing singles of that type in microgravity, 
although it hasn’t been proved otherwise. Today, what is more 
promising, in the opinion of Yu. A. Osipyan, is biotechnology, 
about which V. S. Avduyevskiy spoke. The growth of protein 
crystals in which the macromolecules are very gravity-sensitive 
can be much more efficient in microgravity than the growth of 
semiconductor single crystals, which also show defects in the 
macro- and microirregularity concentration in weightlessness. The 
pharmaceutical industry, which has enormous finances, needs 
samples of perfect macromolecular crystals in order to study the 
mechanism of action of drugs. That is exactly why the production 
of such samples in microgravity could bring quite a return. 


Also of interest are new areas in technologies of metallurgical 
production in connection with the production of eutectic 
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compositions—which, as we know, are very gravity-sensitive—and 
multiphase (dendritic) structure growth. Here, as in the cases 
mentioned above, it is essential to identify active mechanisms 
and assess the influence of the force of gravity. 


The comments of Russian Academy of Sciences Corresponding Member 
R. F. Ganiyev (Institute of Mechanical Engineering, Russian 
Academy of Sciences) were devoted to a discussion of certain 
features of the effects of vibrations in microgravity. In the 
absence of the force of gravity, vibrations can be particularly 
effective in the creation of composite structures or, conversely, 
in the separation of substances. Such vibrational technologies 
are already being used successfully for many production processes 
on the ground. 


Cand. Tech. Sci. A. I. Feonychev (Moscow Aviation Institute) 
spoke on methods for the theoretical description of the 
vibrational components of microaccelerations in spaceflight. He 
noted in particular that a model developed at the Moscow Aviation 
Institute allows for structural rigidity and makes it possible to 
identify changes in the amplitude and dispersion of oscillations 
in microaccelerations. In his opinion, the "convection sensor" 
measurement of microaccelerations, which was discussed in a 
number of papers, is not promising. 


The comments of V. F. Agarkov (Central Design Bureau, Samara) 
touched on the features of the design of special space platforms 
for microgravity research and their certification. At present, 
most of the work that is being done in that area uses unmanned 
vehicles of the Foton and Bion series, which have a level of g/g, 
< 10° for linear accelerations and g/g, < 10° for vibrational 
microaccelerations. In addition, a new generation of space 
vehicles used for production and biotechnology purposes is being 
developed, and they are to have a considerably lower level of 
microaccelerations (down to g/g, < 5°10’ for linear 
accelerations). Among such space vehicles are the NIKA-T 
(production vehicle) and the NIKA-B (biotechnology vehicle). The 
certification of this new generation of space vehicles in terms 
of microaccelerations involves the precision calculation and 
measurement of microaccelerations and the special design of space 
platforms that use devices to attenuate vibrations. 


Dr. Phys.-Math. Sci. S. G. Cherkasov (Scientific Research 
Institute of Thermal Processes of the Russian Space Agency), in 
his comments, said that, although the overwhelming majority of 
reports and comments heard at the session were devoted to 
mechanics problems that involve production processes for 
obtaining various materials in space, hydromechanics in 
microgravity is, on the whole, a much broader research area and 
has a broader field of application than does space technology. An 
example consists in the hydrodynamics and heat-and-mass transfer 
processes that take place in the fuel tanks of aerospace systems. 
Microgravity has a definite effect on the qualitative and 
quantitative nature of the processes in space-vehicle tanks. And 
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here, for fluid mechanics in microgravity as a research area is a 
very broad and important field of application. 


In recent years, our scientists, in studying questions of the 
hydromechanics of weightlessness, have begun to move away from 
problems associated with in-tank processes; but in certain 
questions, that gap can be easily closed. Such is true primarily 
of processes seen in space technology and in-tank processes that 
are close in their physical essence. In cryogenic fuel tanks, for 
example, as in the growth of crystals, natural convection plays 
an important role. In both areas, although on the basis of 
different approaches, substantial results have been obtained. 
That is why coordination of efforts in that area in terms both of 
ideas and organization can have a big effect. 


Cherkasov noted further that, in the field of in-tank processes, 
there is another set of problems that are virtually unstudied in 
fundamental terms. For the long-term storage of cryogenic fuel 
components in space, for example, special cooling systems that 
compensate for background heat flows into tanks are needs, as is 
a related rise in pressure. Such systems must include heat 
exchangers located in the tanks. If a heat exchanger tube being 
cooled is located in the vapor phase of the fuel, then 
condensation will occur on that tube. On the ground, heat 
exchange via condensation is a very effective mechanism for heat 
transfer because of the drip of the film of condensate from the 
heat transfer surface as a result of gravity. But in spaceflight, 
obviously, the picture is very different. 


At first glance, the tube should develop a thick buildup of 
condensate very quickly, and the process should end with that. 
But there are also capillary forces that could induce certain 
kinds of condensate flows. For example, it can be demonstrated 
that a constant-thickness film on a cylindrical tube will be 
unstable. Moreover, excessive capillary pressure can cause the 
film to drip from the tube at points where the free-suriace 
curvature is great onto connective joints and tank walls where 
the free-surface curvature is smaller. Such capillary-force- 
induced flows can, obviously, have a substantial effect on the 
efficiency of operation of a heat exchanger; but it is difficult 
to estimate the rate of the flows, because such flows have never 
been studied. 


Cherkasov emphasized that that example, as well as others, shows 
that the problems associated with the in-orbit production of 
materials represent nowhere near all the problems of 
hydromechanics in weightlessness and microgravity. There are also 
other areas that are important in terms of application, and 
research has been done in those areas and still is being done. 
For that reason, what is quite needed at this stage is 
coordination of efforts both in terms of information on results 
produced and methods developed and in terms of the selection of 
areas for further research. 


- 31 - 





Prof. V. S. Zemskov (A. A. Baykov Institute of Metallurgy, 
Russian Academy of Sciences) noted in his remarks the importance 
of the issues discussed in the papers read at the science session 
regarding the theoretical and experimental research in the field 
of hydromechanics in terms of solving problems associated with 
the growth of uniform single crystals of semiconductors. The 
processes associated with the growth both of volumetric single 
crystals and of layers of semiconductors are accompanied by 
complex heterogenous processes that take place at the interfaces 
of multicomponent solid, liquid, and gaseous phases. Quite 
naturally, of course, one cannot do without the methods and 
mathematics of hydromechanics in the review and analysis of those 
processes. 


Zemskov reported that the very first in-orbit production 
experiments showed that crystallization processes have a high 
gravity-sensitivity. For example, int the Soyuz—Apollo "Universal 
Furnace" experiment, singles crystals of a solid solution of 
germanium—-silicon alloyed with antimony turned out to be less 
uniform in terms of composition and more ideal in terms of 
structure than were similar crystals produced on the ground. Such 
unexpected results in the production experiments served as an 
incentive to formulate and develop theoretical and experimental 
research in the hydromechanics and processes associated with heat 
and mass transfer in weightlessness. A great deal has been done 
in that area. The principal achievement and advances associated 
with that research were reported and discusses at this science 
session. 


Zemskov also noted that, to date, little attention is being 
devoted to the use of the ideas and mathematics of 
physicochemical hydrodynamics to analyze the processes associated 
with the growth of single crystals and semiconductor layers. The 
uniformity of the composition of the single crystals and layers 
depends not only on the nature of heat and mass transfer and on 
the magnitudes and directions of the vectors of the flows in the 
melt (or gaseous phase), but also on the physicochemical 
parameters of specific components of the crystallizing system, 
e.g., parameters such as the coefficients of distribution of the 
components between phases and the coefficients of diffusion of 
the components in a melt. 


Work done aboard Kosmos satellites and in ground-based 
laboratories by the Russian Academy of Sciences Institute of 
Metallurgy in a study of the processes associated with oriented 
crystallization in ampules containing a germanium melt doped with 
various additives has shown that the uniformity of ingots grown 
depends both on the conditions of mass transfer in the melt and 
on the magnitude of the coefficients of the distribution of the 
components between the melt and the solid phase. For each 
component, there are gravity conditions in which the desired 
uniformity of distribution of components in the ingot can be 
achieved. 
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Zemskov also reported that recent work quite rightly devotes much 
attention to the effects produced by an instability of 
microaccelerations (or vibrations) on the processes associated 
with heat and mass transfer during in-orbit crystallization. 
There are, however, experimental data that indicate that the 
constant microacceleration component (linear acceleration) on a 
space vehicle exerts a very strong influence on the uneven 
distribution of the components in a melt. 


Specifically, research has been performed at the Russian Academy 
of Sciences Institute of Metallurgy on the distribution of 
components in the melts of an indium—bismuth-antimony system in 
which the linear accelerations of the force of gravity are 10° g, 
and 10°° g,, and g, is the acceleration of the force of gravity on 
the ground. Experiments at 10° g, and 10° g, were performed 
aboard Salyut 6—Soyuz when the station was in gravity-stabilized 
mode and was in a special rotation. 


In order to achieve total homogenization of the system under 
study, the linear microacceleration could not exceed 10° g,. 
Analysis of the transfer processes in weightlessness would be 
incomplete without consideration of the linear microacceleration 
component. A complete description of the gravity conditions 
aboard a space vehicle requires measurements not only on the 
vibration components, but also of the linear microacceleration 
component. The data cited were presented at the science session 
in the form of poster reports. 


Dr. Tech. Sci. B. G. Zakharov (branch of the Institute of 
Crystallography, Russian Academy of Sciences) noted in his 
remarks the inadequacy of the link between basic-science research 
and the application of that research to the processes that take 
place in melts during crystallization. It is for that reason, 
specifically, that the requisite basis for building a physical 
model and performing systematic calculations does not exist. That 
is why, in the formulation of a program of space-based materials 
science, work in which studies of the processes that occur on the 
ground and in microgravity are closely coordinated with each 
other should receive special support; and methods for producing, 
say, crystals should be based on standard models and gear, as 
should a comparison of results. The focus of the research needs 
to be shifted to detailed preparation and the conduct of a 
ground-based segment of the experiments and calculations that, as 
a scientific base, would ensure the success of the in-orbit 
experiments in the future. 


It is Zakharov’s opinion that, at present, of great interest for 
materials science and microelectronics is research involving the 
perfection of the technologies used for the production of single 
crystals of gallium arsenide in which the capabilities for 
improving the microuniformity of crystal structure and properties 
are linked to a reduction of the deviations from stoichiometry in 
liquid phase. Of course, that must be reproduced in microgravity; 
and when those experiments are set up, the ground-based 
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processing of the results, with detailed mathematical modeling 
and comparison of all possible alternatives for improving the 
technology (among which, microgravity is only one possibility) 
must not only be mandatory, but must also predominate. The center 
created at the branch of the Institute of Crystallography in 
Kaluga will be doing that. 


Prof. V. I. Polezhayev (Institute of Problems of Mechanics, 
Russian Academy of Sciences) noted in his comments that, right 
now, space-based materials science and mechanics of 
weightlessness represent one of the leading areas in space 
research, and, as a result, national and international space 
agencies abroad have divisions that coordinate that field of 
activity, with the level of financing representing anywhere from 
tens of millions of dollars a year to hundreds of millions of 
dollars. The payload market in the United States, Japan, and 
Western Europe is quite large and is growing, particularly 
because of countries that are striving to develop high-level, 
science-intensive technologies (such as Taiwan, South Korea, 
Turkey, and China). 


In Russia right now, there are no customers wealthy enough for 
such operations, because domestic technology doesn’t yield high 
profits, and finances for basic-science research are limited. But 
Russia’s aerospace systems, if operations were efficiently 
organized, could compete successfully with foreign systems. But 
the payload is delivered, as a rule, in "black box" form, i.e., 
without the involvement of domestic specialists in the various 
subject fields, which narrows the potential market of users and 
does little to develop the scientific areas that could later 
support that field of activity of the aerospace scientific- 
production associations. 


Polezhayev also noted that, for more than 15 years, the USSR’s 
programs for space-based materials science and the mechanics of 
weightlessness were conducted on a broad front, with centralized 
financing through the head organizations TSNIIMASH (of the 
Ministry of General Machine Building) and the Space Research 
Institute (of the Academy of Sciences), which, as a rule, copied 
the trends associated with the development of microgravity 
science in the West. For example, the "legend" of the incredible 
profits coming from the biotechnology industry in space was 
disseminated and outfitted materially, as a result of which more 
electrophoretic chambers were built in the USSR (where there was 
still no advanced ground-based industry) than in the United 
States, where that area of science was cut back. 


In Russian at present, the use of orbital stations and probes for 
work involving space-based materials science is unavoidably 
linked to access to the international payload market, with all 
the attendant consequences: competition and the need to maintain 
a given technical level, devoting particular attention to the 
relatively distant future (international space laboratories take 
7-10 years to prepare), painstaking appraisals of the leading 
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research (especially analysis and interpretation of results), and 
the selling of technical in-orbit payload capabilities and the 
advertising of such. One of the top priorities is the maintenance 
of controlled conditions of long-duration spaceflight in 
weightlessness, to include conditions in which other factors have 
an effect (the vacuum and solar radiation). It should be borne in 
mind that, under the new conditions, scientists who worked before 
on state-supported programs involving space-based materials 
science and then became partially or completely unemployed have 
the alternative of entering the international science-product 
market. 


For that reason, research must be aimed at developing 
recommendations for the design and application of orbital 
laboratories for space-based materials science, at appraising the 
work that has been done, at fleshing out areas involving the most 
effective use of microgravity conditions on orbital complexes, at 
conducting a search for expert information on the most efficient 
commercial use of the conditions of microgravity, and at 
preparing and performing experiments involving microgravity 
during joint international missions. 


As for the areas of research themselves, it would be wiser for 
the near term to have a research program that is geared to 
efficient payload use, with the application of a quantitative 
body of mathematical and physical modeling for analysis and 
interpretation of results, and that is narrower than the broad 
programs of space-based materials science and technology that 
used to have budget financing. 


Expert analysis and interpretation of experiments performed for 
the production of materials and substances aboard orbital 
stations should take into consideration actual flight conditions, 
should be based on already developed techniques and program 
complexes, and should contain an analysis of gravity sensitivity, 
as well as conclusions and recommendations on precluding defects 
in future experiments. Only via that route can conclusions on the 
prospects of industrial and semiindustrial production be made. It 
is extremely important to organize international symposia and 
take part in the work of commissions on microgravity for the 
purpose of attracting payload customers from among foreign 
science centers for microgravity, space agencies, and space 
companies. 


In summing up the round table discussion, Academician G. G. 
Chernyy said that this science session displayed the state of the 
work being done with regard to microgravity mechanics and a 
number of its applications, and in that regard, the goal that was 
set before the session was achieved. Obviously, the field of 
applications of microgravity mechanics is extremely broad and is 
not limited to the in-orbit crystal-growth that was referred to 
repeatedly in this session. 


Of extreme interest are problems associated with condensation in 
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weightlessness as they relate to the possibility of using phase 
transition to create efficient heat-exchangers (or heating 
tubes). In weightlessness, because of the absence of the pull of 
the Earth’s gravity, one can use, for example, electrical fields 
to create centrosymmetric mass-force fields to model the 
atmospheres of the planets and stars. 


Immense problems in mechanics occur when one examines the 
processes that take place at the crystallization front, which the 
phenomenological models that are widely used cannot fully 
describe. That is also true of dendritic crystallization. The 
applications of physicochemical hydromechanics in that field are 
poorly developed by comparison with problems involving, for 
example, supersonic gas dynamics, where large groups of 
specialists have worked in that field. For that reason, the 
science potential assembled in mechanics can be applied 
successfully in such a comparatively new field of mechanics as 
microgravity hydromechanics. 


In conclusion, Academician G. G. Chernyy remarked the extremely 
intense development of a whole complex of microgravity sciences 
abroad, and he included many of the microgravity hydromechanics 
issues covered in this session. Foreign companies are displaying 
a great deal of interest in performing their experiments aboard 
our manned and unmanned vehicles-the Foton satellites in 
particular-—proposing black-box payloads. But there are also 
proposals for joint research: upon visiting the Russian Academy 
of Sciences Department of Problems of Machine Building, 
Mechanics, and Control Processes, for example, H. Walter, the 
chief research associate for the microgravity division of the 
European Space Agency, made a number of specific such 
suggestions. 
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Seminar on Numerical Methods Used in Problems Involving Heat-and- 
Mass Transfer, led by V. I. Polezhayev, L. A. Chudov, and G. S. 
Glushko 


957Q0020E Moscow IZVESTIYA ROSSIYSKOY AKADEMII NAUK MEKHANIKA 
ZHIDKOSTI I GAZA in Russian No 5, Sep-Oct 94 pp 171-175 


[Article consists of abstracts of papers given at seminar, under 
the rubric "Institute of Problems of Mechanics, Russian Academy 
of Sciences, Seminars"; UDC 531/534:0.61.3] 


[FBIS text] 


26 Oct 1993. A. S. Demidov (MGU, Moscow). Mathematical experiment 
involving a study of heat-and-mass transfer in the zone of 
evaporation of a heating tube that operates on the basis of the 
principle of an inverted meniscus. The paper presents a 
theoretical study of heat-and-mass transfer when there is 
evaporation in a moistened capillary structure that comes in 
contact with a wall that is heating up. The wall has grooves 
along which the vapor is removed. The set up and results of the 
numerical study of the mathematical model are given: the model is 
a stationary two-dimensional boundary-value problem with a free 
boundary that corresponds to the surface of the phase interface 
within the capillary structure. The generation and growth of the 
vapor zones (bubbles) are tracked. 


15 Nov 1993. A. M. Frank (Computer Center, Siberian Department, 
Russian Academy of Sciences, Krasnoyarsk). Discrete models of an 
incompressible fluid. A number of discrete models of a nonviscous 
incompressible fluid based on the Hamilton principle and using 
regular LaGrangian systems, Dirichlet systems, and an arbitrary 
movable multiangle system, as well as a model of free particles 
based on the Gauss principle (PMTF, 1989, No. 3; ZHURN. VYCHISL. 
MATEMATIKI I MAT. FIZIKI, 1988, No. 11; PMTF, 1993, No. 5). The 
results are given for calculations involving problems of wave 
hydrodynamics: generation of consolidated waves, transformation 
of such waves at depth discontinuity, run-up against wall and 
gently sloping shore, interaction of waves with obstacles, Mach 
reflection of consolidated waves off oblique wall, collapse of 
surface waves. 
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22 Nov 1993. V. A. Sudarikov (Science Center for Neural 
Computers, Russian Academy of Sciences, Moscow Aviation 
institute, Moscow). Neural-network algorithms for solving 
problems of mathematical physics. Cellular neural networks for 
modeling physical systems described by equations in partial 
derivatives are examined. Neural-network algorithms are built for 
solving elliptical and parabolic equations in partial 
derivatives. For a Poisson equation in a two-dimensional space, a 
method is formulated that uses a recurrent, modified 10-digit 
system that uses various activation functions. For an equation of 
thermal conductivity, a neural network is built that makes it 
possible to solve the equation simultaneously and collaterally. A 
problem is examined in which optimal curvilinear grids are built 
on a self-teaching neural network. 


29 Nov 1993. G. S. Glushko (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Expressions for components 
of tensors of coefficients of turbulent viscosity and thermal 
conductivity in two-dimensional flows. The use of an advanced 
mathematical model of turbulent transport processes to derive 
expressions for the components of the fourth-dimension tensor of 
the coefficient of viscosity and the second-dimension tensor of 
the coefficient of thermal conductivity is examined. The 
expressions obtained are demonstrated, as is the type of system 
of Navier-Stokes differential equations for two-dimensional (on 
average) turbulent flows. 


A system of two-dimensional equations obtained via averaging for 
one of the spatial directions is derived: equations for velocity 
components, and a vorticity equation. The additional members that 
appeared in the vorticity equation are discussed, and a problem 
involving the determination of the coefficients of vorticity 
transport is formulated. The prospects of the further development 
of the mathematical model of the turbulent transport processes 
are discussed, as are possible applications. 


6 Dec 1993. A. A. Gorbunov (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Numerical modeling of the 
hydrodynamics of a self-gravitating volume of relativistic gas. 
The spatial motions of a relativistic gas and the evolution of 
the thermodynamic parameters in a centrosymmetric force field 
that arise when adiabatic equilibrium in some vicinity of the 
central point is disturbed are studied numerically. The defining 
equations are written in a Cartesian system of coordinates. The 
initial and boundary conditions of the problem are discussed. A 
finite-difference numerical method is elaborated. The stability 
of the difference systems used is studied on a qualitative basis. 
The possibility of using conveyer and multiprocessor computations 
is examined. 


13 Dec 1993. D. S. Pavlovskiy (Institute of Problems of 
Mechanics, Russian Academy of Sciences, Moscow). Analysis of the 
stability of a convective flow. A problem is examined that 
involves convection in elongated flat layers with a free surface 
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on which temperature-dependent surface tension is active. The 
convection comes about as a result of a longitudinal temperature 
gradient along the layer. A problem involving the stability of a 
plane-parallel flow is worked through. Methods of slightly 
nonlinear analysis are used to find secondary flows when there 
are weak supercriticalities. 


20 Dec 1993. V. A. Frost (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Equation for probability 
distribution function of concentration, and problem of modeling 
turbulent streams with variable density. An equation for the 
density function of concentration and an equation for the joint 
density function of velocity and concentration are examined as 
the basis for constructing a mathematical model of. turbulent 
combustion. Primary attention is focused on a description of 
turbulent variable-density streams in the presence of volumetric 
forces (gravity, mean pressure gradient). The features of the 
behavior of turbulence in a pulse combustor are examined for 
purposes of illustration. 


20 Dec 1993. S. T. Surzhikov (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Mathematical modeling of 
the subsonic motions of a radiating gas. A method based on the 
solution of vortex-—current function-variable Navier-Stokes 
equations, with the introduction of a new function that ensures 
the existence of a nonstationary analogue of the current 
function, is examined for the calculation of the nonstationary 
subsonic motions of a viscous, compressible, heat-conducting gas. 
Problems are solved that involve the appearance of vortex 
structures in optical plasmatron chambers, the motion of a 
radiating gas in a subsonic nozzle, and the dynamics of radiating 
burn volumes in a field of gravity. 


24 Dec 1993. I. A. Kryukov (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Program package for solving 
two-dimensional nonstationary Eulerian equations. The theoretical 
bases for numerical methods used for calculating nonviscous, 
compressible flows are examined: TVD and ENO retrieval 
algorithms, methods for solving problems involving the breakdown 
of arbitrary discontinuity, methods for sampling equations for 
the flight coordinate, etc. Examples are given of equations of 
nonstationary problems of gas dynamics in one- and two- 
dimensional conditions. 


27 Dec 1993. A. I. Prostomolotov (Institute of Problems of 
Mechanics, Russian Academy of Sciences, Moscow). Chokhralskiy 
study of convective flows due to low-frequency crystal 
vibrations. Experimentally studied is the evolution of the 
spatial structure of a current with growth of the vibrational 
Reynolds ..mber and varying height of the subcrystal fluid 
column. Experiments are conducted on model fluids in a broad 
range of Re numbers. Detected in cylindrical crucibles of regular 
structure were surface currents and intense large-scale flows in 
the volume of the fluid as a result of low-frequency vibrations 
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imposed on a solid coming in contact with the fluid. Analysis 
shows that a controlled vibrational action directed at the 
surface of the fluid can reduce considerably the uncontrolled 
effect of thermal capillary or concentration-capillary convection 
on the transport processes when crystals are grown with the 
Chokhralskiy method. 


28 Dec 1993. N.: A. Verezub (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Use of magnetic systems in 
the M technology for growing single crystals of silicon from 4 
melt. The optimality of the hydrodynamic picture of a current, 
determined for a low level of oxygen in single crystals of 
silicon by the requirement that the subcrystal circulation of a 
melt be separated from the crucible walls, is chosen as the basic 
criterion for optimizing the parameters of the magnetic system 
(configuration and intensity of magnetic flux). A cycle of 
parameter studies is performed for the purpose of analyzing the 
effects that rotating and permanent-magnet fields of varying 
configuration and intensity have on the flow of a silicon melt in 
the Chokhralskiy method. Hydrodynamic processes are modeled on 
the basis of the solution of three-dimensional equations of 
hydrodynamics and heat-and-mass transfer via the finite- 
differences method. At the same time, a study is performed of the 
growth of single crystals of silicon in magnetic fields under 
factory conditions. Satisfactory agreement with experimental data 
is noted. Specifically, in one of the experiments, the rate of 
rotation of the melt, measured on the basis of the rate of 
rotation of the residues of charge, is 3.4 cm/s; while the rate 
calculated from the experiment parameters is 4.65 cm/s. 


11 Jan 1994. M. K. Yermakov, S. A. Nikitin, D. S. Pavlovskiy, V- 
I. Polezhayev, S. L. Sternin (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Problems of mechanics 
involving certification of space vehicles designed for space- 
based materials science and biotechnology. New problems of 
mechanics are examined that involve the entry of domestic space 
hardware onto the international payload market in the field of 
space-based materials science and biotechnology. Results are 
provided of the analysis of low-frequency and high-frequency 
microaccelerations and their effect on convective processes in 
the growth of crystals and the separation of substances. Examples 
are given of the calculation of heat exchange in space-vehicle 
compartments during orbital flight. 


14 Jan 1994. K. G. Dobovik (Institute of Problems of Mechanics, 
Russian Academy of Sciences, Moscow). Appearance of oscillations 
in oriented crystallization in weightlessness. One possible 
mechanism of the appearance of hydrodynamic oscillations due to 
double Maragony diffusion (interaction of thermal capillary and 
concentration-capillary effects) is examined. The possibility is 
demonstrated for the existence of three different flow regimes. 
Two types of oscillations are identified. A means of problem 
simplification is proposed. 
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31 Jan 1994. A. I. Fedyushkin (Institute of Problems of 
Mechanics, Russian Academy of Sciences, Moscow). On certain 
problems of convective mass-and-heat transfer. On the basis of 
Navier-Stokes, the set up and results of the numerical solution 
of the following problems are examined: convective heat-and-mass 
transfer in the crystallization of eutectic alloys 
Si/Al—-Ni/Al-Cu/Ge, the mathematical modeling of convective heat- 
and-mass transfer in nonrectangular regions, and the mathematical 
modeling of convective heat-and-mass transfer in various 
gravitational fields. 


7 Feb 1994. S. F. Timashev (Karpov Scientific Research Chemical 
and Pharmaceutical Institute). Principles of Flicker-Noise 
Spectroscopy. It is shown that the well-known universality of 
flicker-noise (1/f-noise) that is seen in various physicochemical 
and natural processes is the result of effects of alternation in 
the evolution of complex dynamic systems. General expressions are 
produced for the spectral density of the autocorrelator of the 
dynamic variable in such evolution, which enables the use of the 
parameters of 1/f-noise as a source of information on dynamic 
system states and on the change in those states in the course of 
the evolution (flicker-noises, FN spectroscopy). The possibility 
is discussed of using FN spectroscopy in various regions: 
electronics, materials science, geophysics, astrophysics, 
medicine, ecology, genetics, and economics. 


14 Feb 1994. V. D. Golyshev (All-Russian Scientific Research 
Institute of Synthetic and Cleaning Agents). Experience in the 
development and use of methods of growing large crystals in a 
one-dimensional thermal field in a melt. The thermophysical 
formulation of a problem for the development of a growth method 
whose thermal conditions do not depend on the transverse cross- 
section of a growing crystal is described. Descriptions are given 
for designs of two arrangements for two different types of 
methods. Six versions of the method are analyzed, as are the 
results of their experimental testing. Data are cited for 
experiments involving the precision of determination of thermal 
fields in the growth of crystals, and results are given for the 
analysis of the link between mica defects and thermal conditions. 


14 Mar 1994. N. V. Nikitin (Institute of Applied Mechanics, 
Russian Academy of Sciences, Moscow). Numerical study of 
turbulence and transport in circular-section tubes. A problem of 
numerical modeling of turbulence and transport in circular- 
section tubes is studied. Three-dimensional Navier-Stokes 
equations are solved in a cylindrical system of coorainates (x, 
r, 6). Two formulations of the problem are discussed: one with 
conditions of periodicity in the direction of the main flow, and 
one without such conditions. The results obtained in both 
formulations are reported. Method of approximation: finite- 
differences in a radial direction; spectral (Fourier series); or 
finite-differences in the x and 6 directions. Time integration is 
done with the semi-implicit Runge-Kutta method with automatic 
step choice. 
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Periodic boundary conditions are set up in the modeling of 
advanced turbulence in the direction of the main flow. More than 
40 calculations are performed at Reynolds numbers of up to Re = 
10* with a spatial of from 32 x 32 x 32 to 64 x 64 x 64. The 
effects of algorithm parameters on computation results are 
studied. Calculations were continued until statistically 
stationary regimes of motion were established and then until 
reliable results of averaging were obtained. A detailed 
comparison was made with well-known experimental data, and the 
comparison demonstrated a good overall agreement between the 
calculated currents and the actual turbulent movements in long 
tubes. For modeling laminar-turbulent transport, outlet boundary 
conditions were proposed that made it possible to complete a 
problem on the spatial development of turbulence. Calculations 
were performed that modeled the transport in tubes of finite 
length: x,,,/Yr. < 120 (x,,, is the length, r, is the radius of the 
tube). In those calculations, a spatial resolution of 512 x 32 x 
32 (in the directions x and r, respectively) was used. The 
results obtained qualitatively and quantitatively agree with the 
experimental observations. 


21 Mar 1994. P. N. Babishchevich (Institute of Mathematical 
Modeling, Russian Academy of Sciences, Moscow). Difference 
methods for solving problems of hydrodynamics and heat-and-mass 
transfer in current function—vorticity variables. Difference 
methods for solving two- and three-dimensional problems involving 
the dynamics of a viscous incompressible fluid and convection in 
complex design fields in current function-—vorticity variables on 
nonuniform rectangular grids are examined. In irregular fields, 
computational algorithms for continuous counting are built on the 
basis of the method of hypothetical fields-solution of the 
problem in an expanded regular field. The main features of the 
difference systems that were built are as follows: the difference 
systems were linearized; at every time layer, self-coupled and 
nonself-coupled elliptical grid problems were solved. The 
algorithms for realizing the boundary conditions for the 
vorticity were noniterative and absolutely stable. In special 
approximations of the convective members and in the 
disintegrations, the linearized systems were absolutely stable: 
the appropriate a priori estimates were made for the difference 
solution for all grids. 


Iterative solution of self-coupled and nonself-coupled elliptical 
grid problems was performed for realizing the systems. 


- END - 
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